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Abstract 



Universal scaling in the power-law size distribution of pelagic fish schools is estab- 
lished. The power-law exponent of size distributions is extracted through the data 
collapse. The distribution depends on the school size only through the ratio of the 
size to the expected size of the schools an arbitrary individual engages in. This 
expected size is linear in the ratio of the spatial population density of fish to the 
breakup rate of school. By means of extensive numerical simulations, it is verified 
that the law is completely independent of the dimension of the space in which the 
fish move. Besides the scaling analysis on school size distributions, the integrity of 
schools over extended periods of time is discussed. 
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1 Introduction 



Pelagic fish commonly cruise as a school. As they migrate about in a limited 
space, interaction between schools can occur so that two of them encounter and 
aggregate, or a large school splits itself into two smaller schools, or more. The 
continuation of these interaction processes of f ission and fusion may eventu- 
ally lead to a fat-tailed school-size distribution. iBonabeau and Dagor 

311991) 



found power law in school-size distributions of tropical tuna. T he pow er-law 
distributions of school sizes are quite general for pelagic fishes flNiwal .E998l. 



and quantitative analyses are now in progress (Bonabeau et al., 1998, 1999; 
Niwa, 2003). 
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A model of school aggregation proposed by Bonabeau and Daeorn ( 19951 ) is 



based on a physical model of partic le aggregation, i .e. open systems showing 
power-law cluster-size distributions ( Takavasu . 19891 ). Bonabeau and Dagorn's 



(1995) model is briefly sketched as follows. Simulated fish schools move be- 
tween sites on coarse-grained zones of space, and aggregate when they meet. 
School splitting is replaced by (re)injection process: a certain fraction of each 
school is separated from the school; the fish that have left their school are rein- 
jected as individuals (i.e. one-sized schools') into the systena, wh ile the total 
number of individuals is kept fixed. iBonabeau and Dagorn obtained a 



truncated power law with exponent —3/2 in the mean-field case of the model, 
in which schools can move from any site to any other site. 

The mean-field assumption might not always be adequate, because there are 
spatial constraints on movement. Although the ocean is three-dimensional, 
fish may not fully use their spatial environment. Pelagic fish movement gener- 
ally takes place in a horizontal (two-dimensional) space. They are limited in 
depth by physiological constraints, and do not dive into the deeps. Addition- 
ally, pelagic fish schools are concentrated in the vicinity of t he front wh ich is 
the contact zone and collision line of two oceanic currents (Uda, 1938[ ). For 



instance, skipjack schools, Euthynnus vagans (Lesson), aggregate in t he in- 



terfacial region between the cold subarctic and warm subtropic waters (jUda . 



1936| ). They are constrained to move effectively in a one-dimensional space. 



In the presence of a fish aggregating device (FAD), which is a drifting log or 
an artificial device designed to attract fish, pelagic fish do not make full use of 
the three-dimensional oceanic space and are concentrated in the vicinity of a 
FAD (a point), so that the space dimension effectively decreases to less than 
one. 

Bonabeau et al. (1998, 1999) predicted spatial effects, i.e. the dimensional re- 
duction that the power-law exponent of fish school-size distributions is modi- 
fied from the mean-field case of their model when the dimension of the space 
is taken into account: the absolute value of the exponent, in a version of their 
model on a d-dimensional lattice (schools hop to neighboring sites only), de- 
creases when the space dimension d decreases. 



Recently, Niwa ( 2003| ) analyzed some existing data for various species in terms 



of a school-size histogram of the population of fish. fish swimming together 
form an A^-sized school. The school-size distribution W{N) is proportional to 
the observed number of A^-sized schools; the size histogram of the population 
is then represented by P{N) = NW{N), which is proportional to the fraction 
of fish in A^-sized schools to total population. He reported that the distribu- 
tion P is an exponentially thin-tailed distribution, and, the distribution W 
follows a power-law decay with exponent —1 and is truncated at a cut-off size. 
A simple stochastic-differential-equation model was proposed to explain the 
observed power-law behavior, and the predictions of the model were found 
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to be consistent with empirical data. A remarkable feature is "scaling": all 
the empirical distributions collapse onto a single curve if the data are plot- 
ted in terms of scaled coordinates with the mean value of a histogram P{N), 
for all various species and for all environmental factors including the above 
mentioned space dimension. Note that the power-law distribution W does not 
have a well-defined mean; contrary the rapidly decreasing distribution P has 
a well-defined mean. 



The empirically determined power-law exponents of school-size distributions 
for specific data sets range from 0.7 to 1.8 (Niwa, 1998; Bonabeau et al., 
1999). Fat-tailed school-size distributions are necessarily truncated because 
the population is finite. This truncation of the power-law regime might lead 
to a "wrong" estimation of the exponent. I make use of the data collapse to 
extract the "right" exponent. I here propose a finite-size scaling (FSS) form 
( Binder and Heermann . 1988[ ) for the size distribution W{N) on the assump- 



tion that the distribution decays with the truncated power-law form with 
biologically universal exponent. The power-law exponent of distribution W 
is determined by experimentally fitting the FSS relation to achieve the best 
data collapse. It is investigated by simulating the school system of pelagic fish 
whether or not the dimension of the space in which fish swim is relevant to the 
power-law exponent. FSS for fish school sizes is elucidated through the com- 
petition between two processes in the interacting school system: aggregation 
and splitting of schools. 

In addition, it is numerically investigated how long a school stays together, 
i.e. neither merges with any other schools nor breaks up. The behavioral algo- 
rithr as governing school f o rmati on and dynamics have been extensively studied 
[e.g. Inada and Kawachi ( 20021 ) and references cited therein ] . Stu dies of the 



i ntegr ity of schools, however, are very f ew [e g. iLester et all (|l98,^ or iBavhfj 



( 1988[ ) for experimental studies; Niwal ( 1996f ) for a modeling approach]. Nu- 



merical simulations of fish-school aggregation suggest conjectures about real 
situations that could be tested by observations. 



2 The Data 



The enlarged sets of the data in iNiwal (|2003[ ) for school-size distributions of 
pelagic fishes are analyzed (summarized in Table 1). The ways of estimating 
school sizes of pelagic fishes were catch per set by a purse seine or acoustic 
surveys. Catch-per-set data are expressed in school weight (in metric tons). 
Acoustic-survey data are expressed in dimensional size of a school (e.g. verti- 
cal thickness in meters), which can be reduced to the biomass in a school: the 
school biomass is proportional to the vertical cross-section, the square of ver- 
tical thickness, or the square of diameter of a school (Squire, 1978; Anderson, 
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Table 1 

Species analyzed. 





Species 


cut-off size*^ 


Data sources 


A 


Northern anchovy 


15.67 • AN 


Smith fl970)^ 




Engraulis mordax ^ 




Acoustic survey 


• 


Japanese sardine 


6.45 • (AiV/2) to 

V If 


Kara, (1990). 




Sardinops melanosticta 


IZ.ll ■ {AN/2) 


22 acoustic surveys 


□ 


Tropical tuna'^ 


11.69 • AA^ 


Bonabeau et al. fl999). 




(free swimming) 




Data from fisheries 





Tropical tuna'^ 


4.80 • AN 


Bonabeau et al. fl999'). 




(caught in the vicinity of FADs) 




Data from fisheries 


o 


Herring Clupea harengus 


7.24 • AN to 


Eeid et al. f2000l. 






10.88 • AA^ 


4 acoustic surveys 



^The cut-off size of power-law distribution is calculated by Eq.(9). A A denotes 
the class width of frequency data of school sizes, and will be omitted, 
''possibly including Trachurus symmetricus, Sarda chiliensis, Scomber japonicus, 
and Sardinops sagax. 
'^Data are cited in Anderson fl981 1. 

"^Three species ( Thunnus albacares, Katsuwonus pelamis, and Thunnus obesus) 
are mixed. 

1981; Misund, 1993; Niwa, 1995; Misund and Coetzee, 2000). 

Let R be the dimensional size (vertical thickness or diameter acoustically 
measured) of a steady moving A^-sized school. The following relation between 
dimensional and social sizes holds in a statistical sense: 

R = (constant) x N'', (1) 

where u = 0.5. The prefactor is supposed to be constant for each data set (i.e. 
each survey). It may depend on the species and vary with regions, seasons 
and years in which fish schools are surveyed. The acoustic-survey data are 
transformed into a social size histogram as follows 

W{N)dN oc W{R)R^^''^^dR, (2) 

where W denotes the school-size distribution density. 

The data are given by the set { (Aj, Wi) \i = 1, 2, . . . , n}. WiAN reads the fre- 
quency of school sizes which lie within the i-th class [Aj — AN/2, Ni + AN/2), 
where AA^ denotes the class width and the i-th class mark is given by Aj = 
{i — 0.5) • AA^. From now on, to simplify the expression, we will omit AA^ (or 
AN/2 for Japanese sardine) in mathematical formulae for processing empiri- 
cal data, that is, such unit of school size is introduced as AA^ = 1 (AA^ = 2 
for Japanese sardine). A school of unit size contains a certain number of in- 
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Fig. 1. Comparison of the empirical data with Gaussian distribu tion (broken line). 
Japanese sardine schools were surveyed in the summer in 1982 dHaral . llQQOl 'l. The 
trial Gaussian distribution has the variance M2 equal to the second moment of em- 
pirical distribution (M2 = Yji^i^i/Hi^i = 46-8). The normalization of Eq.(17), 
Yl^=i-^i^i — {^)p-> is chosen, where {N)p is given by Eq.(9). The Gaussian dis- 
tribution is consistently normalized. The probability to find sardine schools ranging 
from 18 to 20 meters in vertical thickness (denoted by pig) is 2.0 x 10~^ from the 
data. The Gaussian distribution gives pig = 8.1 x 10"^'^. The social size of such 
sardine schools as € [18,20) in meters are tra nsformed int o A'^ G [42,50). The 
solid line shows the size distribution proposed by Niwal ( 20031 ). with {N)p = 11.7 
calculated from the data. Plotted on semilogarithmic scale. 



dividuals, for instance, AA^ should be a number of fish correspond to 1000 
kilograms for tropical tuna. 

As with fisheries such as fishing by a purse seiner, such major events (i.e. 
large catches) can be regarded as unusual events of great magnitude lying on 
the tail of a distribution comprising events that are mostly of much smaller 
magnitude. Figure 1 shows the school-size distribution for Japanese sardine 
Sardinops melanosticta, which was from the acoustic survey off southeastern 
Hokkaido covering the peri od July 30 - August 6 in 1982 ||Hara ffl99Qi) : the 



same data are available in iHaral (|l984^ ]. A traditional, widely used Gauss 



statistics says that finding sardine schools ranging from 18 to 20 meters in 
vertical thickness should only occur about once every 10^ detections of schools. 
In other words, it is not the real world! Aquatic observations actually say that 
finding such schools occurs about once every 500 detections. The probability 
that such schools are found is 10^ times large! 
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2.1 Finite-size Scaling 



The size distributions of pelagic fish schools follow a power law W{N) oc 
up to a cut-off size 

L= ^'=' \ \ (3) 

In order to characterize the effects of the finite population size on the trun- 
cation of power-law distribution, a finite-size scaling hypothesis is used: the 
distribution depends on N only through the ratio N/L"^, 

W{N; L)dN = L-^F (n/L^) d (n/L^) , (4) 

where F is a universal function independent of system (population) size, and 
Wi = Wi/ Er=i N^Wi. The pref actor L'^ is required to ensure the normaliza- 
tion 

n 

i=l 

School sizes are assumed to obey FSS with biologically universal exponents A 
and B for a wide spectrum of both pelagic species and environmental condi- 
tions. The normalization Eq.(5) together with the postulated universal func- 
tion F{x) gives 

L^^-^x^F(x)dx = 1. (6) 



r 

Jo 



It then follows that 



" = 1 P' 

because all powers of L must cancel out. Prom the FSS hypothesis, it is ex- 
pected that when W{N)L^~^^ is plotted against N/L"^ with correct parameters 
A and B all the empirical data should collapse onto a single curve. The power- 
law exponent of fish school-size distributions, j3, is then evaluated through FSS 
analysis. Besides Eq.(7), if FSS is valid the value of A is 1. 

Let us search for the values of A and B that do the best job of placing all 
the data points on a single curve. To do this, the x-axis is divided into bins 
(Fig. 2a), and the parameters are estimated at values that minimize the mean 
of two-dimensional variance 

e = {ajxf + {ay/yf (8) 

(Lillo et al., 2002, 2003), where a denotes the standard deviation, and x and 
y denote the mean, and the x-axis is chosen as x — N/L^, and the y-axis 
represents y = W{N)L^~^^ . The mean of two-dimensional variance, e, is a 
measure to determine the goodness of collapse. Figure 2b displays the set of 
pairs {A, B) in which the minimum of the mean of two-dimensional variance, 
Cmin, is guaranteed to lie. A good data collapse can be obtained by using 
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Fig. 2. FSS analysis of the empirical data (summarized in Table 1). (a) FSS plot 
of the school-size distribution on double-logarithmic scale. Here y = WL^^^ is 
plotted versus x = N/L^ with A = B = 1. The bins are chosen equally spaced 
on a logarithmic scale as x G ^iQ~i+{k~0-5)/3 ^'^Q-i+(k+o.5)/3-^ ^^^y^ k = 0,1,...,5. 

For each bin two-dimensional variance e is calculated (the rectangle in gray reads 
the interval y ± ay] the mean y is ind i cated by the slit). The solid line is a pre- 
diction of the mean- field theory ( Niwal . 2003l 'l . (b) The region of the yli?-plane in 
which the minimum of the mean of two-dimensional variance exists. The mean of 
two-dimensional variance, e, takes a minimum 6niin for the right choice of (^, .Z?). 
The minimum is found with the precision i.e. width of the minimum, Ae = 10"'^ 
in black region (Ae/emm ~ 2.58 x 10~^). The values of the parameters lie in the 
intervals A = 1.001 ± 0.004 and B = 0.997 ± 0.011. The value of B/A is the esti- 
mate of the power-law exponent /3, and therefore (3 = 0.996 it 0.015. Experimentally 
fitting the parameters A and B to achieve a good data collapse yields the values 
A = B = P = 1, which are precisely the values of the mean-field case of the proposed 
model. « 



the values A ^ 1 and B ^ 1. The power-law exponent derived from the 
FSS collapse is (3 ~ 1. The resulting plot is shown in Fig. 2a. The figure 
confirms the FSS hypothesis, since all the data collapse onto a single curve. 
The school-size distribution follows a power-law decay with exponent —1, and 
is truncated at the cut-off size that equals the mean of another distribution of 
the school system, Pi = NiWi, i.e. the size histogram of the fish population. 
The distribution P decays rapidly for large N, and has a well-defined mean 



Two competing processes, mixing and splitting of schools, determine the school- 
size distribution. The spatial density of fish population conditions the mixing 
rate of schools. The population density of a given species differs among re- 
gions, seasons and years. The breakup rate of schools of a given species may 
vary depending on environm ental conditions, e.g. the lack of food may reduce 
school stability (jMorgani 1198 8). Figure 2a shows that the power-law exponent 
"—1" is robust, while the cut-off size of linear power law varies (Table 1). The 
data-collapse implies that the breakup rate and the population density do not 
affect the power-law exponent but the cut-off size. This indicates that the 
power-law exponent is universal, and hence the overall shape of the distribu- 
tion may result from a simple underlying aggregation mechanism. Moreover, 
as a consequence of the quite good data-collapse, the empirical data do not 
support the dimensional reduction, because the "effective" dimension is re- 
lated to biological or environmental conditions. In the following two sections 
we will numerically examine whether or not space dimension influences the 
power-law exponent. 



3 The Model 



A stationary equilibrium system of a flxed population size (number of indi- 
viduals, denoted by $) is considered in which flsh schools break up and merge 
with other schools. Let us investigate the aggregation process in a discretized 
space and time. There are s sites on the lattice space Q. On every site there 
is at most one school of simulated flsh. Let N{j, t) be the size of the school 
on site j at the t-th time step. At each discrete time step, each school hops 
to a new site or breaks up into pairs of schools of various sizes possible. If 
more than two schools happen to hop onto one site, they coalesce into a sin- 
gle school with the size equal to the sum of the sizes of the incident schools. 
Assume binary splitting independent of school size. Each school with a size 
greater than or equal to 2 splits into a pair of schools with a probability p 
at each time step. The probability p for a school to split per time step (i.e. 
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breakup rate) is independent of its size, and the sizes of splitting schools are 
uniformly distributed: a probability for an A^-sized school to split into M- and 
[N — M)-sized schools is represented by 

K^{N\M, N-M)^ K^{N) = (10) 

for > 2. 

The aggregation process can be represented by the following stochastic equa- 
tion for N{j,t): 

Nij,t + l)^J2Mjkmik,t), (11) 

where Mjkit) is a stochastic variable. If the school on site k does not break, 
Mjk{t) is equal to 1 when the school jumps to site j and equal to otherwise 
{Mik{t) = for i 7^ j). If an A^-sized school on site k breaks up into two 
schools of sizes mN and (1 — m)iV jumping to sites j and j' , respectively, then 
we have Mjk{t) = m, Mjik{t) = 1 — m, and the others vanish {Mik{t) = for 
i ^ where m is also a stochastic variable uniformly distributed in (0, 1) 

on condition that mN is an integer. Mj^it) must be normalized as 

E^.-fcW = l' (12) 

for V/c and Vt, which guarantees the conservation of population. 

In the mean-field case of the model, at each time step, all schools move to- 
wards a randomly selected site, which corresponds to the migration of the high 
potential speed of fish, e.g. free-swimming tuna. They may move to any site 
with equal probability 1/s. In spatial model of school aggregation, schools hop 
to neighboring sites only. In a version of the model on a rf-dimensional lattice, 
they may move to each of 2d neighboring sites with equal probability l/2d. If 
a school on site j breaks, one of splitting schools remains at the site j and the 
other hops to neighboring sites. 

3.1 Finite-size Scaling 

Simulations have been performed with the coarse-grained zones of s = 2^^ 
sites, simulation run = 2^^ time steps, and parameters summarized in Table 2. 
The initial school-system configurations are taken to be random distribution 
of eight-sized schools on the lattice space. Numerical results are shown in the 
next section (in Fig.6b an FSS collapse of W{N) for the two-dimensional case 
is depicted). We search for the values of A and B that place all of the simu- 
lated distributions most accurately on a single curve. The parameters A and 
B derived via the minimization of the measure e to quantify FSS collapse are 
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Table 2 

Parameters used in simulations. Monte-Carlo simulations of fish school ag- 
gregation on the lattice space of s = 2^^ sites with periodic boundary have been 
conducted by using Mersenne Twister ( Matsumoto and Nishimur'al 199^ , a pseudo- 
random number generator, on Scientific Computing System of MAFFIN, Tsukuba 
Japan. 



breakup rate p population <I> (-^)p [ID; 2D; mf] 



♦ 


0.02 


2i4 


1.72; 4.23; 8.97 


■ 


0.02 


2l5 


3.08; 8.57; 17.87 


▲ 


0.02 


2l6 


6.39; 18.57; 38.40 


• 


0.02 


217 


13.44; 39.41; 74.25 


* 


0.02 


2l8 


28.39; 73.83; 157.61 





0.01 


2i7 


19.63; 68.59; 146.61 


□ 


0.03 


2i7 


10.98; 27.00; 48.51 


A 


0.04 


217 


9.13; 20.29; 37.03 





0.1 


217 


5.54; 9.41; 13.94 



^ (N) p was computed from simulation results after 2i7 time steps 
for one- (ID), two-dimensional (2D), and mean-field (mf) cases. 



Table 3 

Summary of exponents A, B, and (3. The exponents are evaluated at values in 
the listed ranges with the precision Ae = 10~^ through FSS analyses of empirical 
data and simulated distributions at the last of run for different population sizes and 
breakup rates in one- and two-dimensional cases. 

A B ^ Ae/e^nin 

data 1.001 ±0.004 0.997 ± 0.011 0.996 ± 0.015 2.58 x 10^3 
ID 0.998 ±0.012 1.020 ±0.026 1.022 ±0.037 2.93 x 10"^ 
2D 0.992 ±0.011 1.019 ±0.026 1.028 ±0.037 1.01 x 10"^ 
mf^ 1 1 1 — 

^ predicted values of the mean- field theory (Niwa, 2003) . 



summarized in Table 3. In the one- and two-dimensional cases, there is a clear 
minimum for A ^ B ^ 1. Therefore, the power- law exponent extracted from 
numerical simulations reads /3 ~ 1, and does not depend on the dimension- 
ality d. This may cause one surprise, because such exponents depend on the 
dimensionality as the critical exponents for scaling behavior in wide varieties 
of physical phenomena, e.g. magnetization, specific heat, size of a polymer, 
and so on. 
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Fig. 3. Plot of (iV)p against the scaled population size (numerical results from 
simulations summarized in Tabic 2). The abscissa represents p/p, the coarse-grained 
population density p (= ^ j s) divided by breakup rate p. Broken lines are fits with 
slopes 0.503 (ID), 1.46 (2D), and 3.03 (mf), respectively, which give the prefactor 
c in Eq.(13). 

4 The Scaling Law 



The cut-off size of power-law distribution (equal to {N) p) results from variable 
individual behavior (i.e. breakup rate p) and fluctuating population density 
(denoted by p). The FSS collapse suggests that the school-size distribution 
Vr(A^) depends on N ^ p, and p only through the variable x = N/ {N) p. Numer- 
ical simulations reveal that (A^) p depends linearly on the spatial population 
density and inversely on the breakup rate (Fig. 3): 

{N)p = cp/p. (13) 

The prefactor c must have dimensions of [length]'^/ [time], and it is indeed 
proportional to the ratio of the coalesce rate to the spatial density of schools, 
as discussed later. 

The frequencies of the amount of A^-sized schools (A^ = 1,2,...,$) at the 
last of simulation run in the two-dimensional case, W{N; p,p), are shown in 
Figs.4a and 5a for a given breakup rate p and a given coarse-grained spatial 
population density p — ^/s, respectively. The frequencies are normalized as 

J2 NW{N) = (14) 

N=l 
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Fig. 4. Simulated school-size distribution (2D) with the breakup rate p = 0.02, 
s = 2^* sites after 2^^ time steps, (a) Frequencies W{N) of A'^-sized schools in the 
two-dimensional simulations with various population sizes summarized in Table 2. 
The straight line illustrates a general trend, i.e. the power-law behavior with expo- 
nent — 1. (b) The frequencies in Fig. 4a re-plotted with W' = W{N)p as a function 
of variable A'^' = Nj p. The scaling causes a shift of the curves in Fig.4a that depend 
on p. All the distribution collapse onto a single curve. Plotted on double-logarithmic 
scale. 



Examining Fig. 4a closely shows that for fixed breakup rate p but increasing 
population density p, the range of the power-law regime increases: the cut-off 
size of power law scales with p as depicted in Fig.4b. The FSS form (4) with 
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Fig. 5. Simulated school-size distribution (2D) with the population density p = 0.5 
(population $ = 2^''), s = 2^* sites after 2^^ time steps, (a) Frequencies W{N) 
of A?^-sized schools in the two-dimensional simulations with various breakup rates 
summarized in Table 2. The straight line illustrates a general trend, i.e. the 
power-law behavior with exponent —1. (b) The frequencies in Fig. 5a re-plotted 
with W' = 1Q~'^W {N)p~'^ as a function of variable A'"' = Np. The scaling causes a 
shift of the curves in Fig.Sa that depend on p. All the distribution collapse onto a 
single curve. Plotted on double-logarithmic scale. 

A — B — 1 together with Eq.(13) reads 

W{N- p,p^ 0.02)dAr = N-^g{N/p)d (N/p) , (15) 
where g{x) is a scaling function. 
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Similar behavior is seen in Fig. 5a when p varies. As p increases, the length of 
the power-law region reduces: the cut-off size of power law scales with 1/p as 
depicted in Fig. 5b, i.e. 

W{N; p = 0.5, p)dA^ = pN-'^g'{Np)d (Np) , (16) 

where g'{x) is a scaling function. Each resulting graph shows that all the 
distributions collapse onto a single curve (Figs. 4b and 5b). 

The same results can be obtained with one-dimensional case of the model, and 
its mean-field case as well. 



4-1 The Unified School-size Law 

Equation (13) implies that the population density must appear as a ratio to 
the breakup rate in the scaling analysis, and the breakup rate vice versa. For 
instance, the following normalization is appropriate for the scaling analysis of 
school based data: 

j2NiWi^{N)p. (17) 

1=1 

The ESS form of Eq.(4) then suggests that the linear dependence of {N)p 
on the population size gives B — 1. By arguments discussed below, the ESS 
form (4) expresses the unified scaling law for school sizes. The scaling behavior 
of the school-size distribution is described by two independent parameters: the 
population density p and the breakup rate p. One of the important results of 
the model is that a single parameter (N) p is sufficient to describe the scaling 
behavior, i.e., the two scaling relations of Eqs.(15) and (16) are unified into 
Eq.(4) because of Eq.(13). 

To understand the unified scaling law for school sizes, it is essential to see what 
determines the cut-off size of power law. The cut-off size is determined by the 
competition between breakup and coalescence of schools. The coalescence is 
regulated by the spatial number- density of schools resulting from a balance 
of the population density and the breakup rate of schools. Equation (13) is 
then fundamental to unify the two scaling laws (15) and (16). Let us trace 
the size change of the school a certain individual (named "A") rides. Let 
be the coarse-grained density of schools (i.e. number of schools per site), 
and, A the probability for a school to coalesce with other schools per time 
step (i.e. coalescence rate). Then the expected increment and decrement of 
the size of "A"-riding school (denoted by Na) are co(p/0)A (accompanied by 
a constant factor cq) and pN\/2, respectively, at each time step of coarse- 
grained simulation [p/cj) reads an (apparent) average size of fish schools]. The 
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expected size-change of "A" -riding school per time step, AA^a, obeys 




ANa = 

The ratio (Xcop/cf))/ {p/2), therefore, gives the expected size of "A"-riding 
school in the stationary state, which is equivalent to the value of {N)p by 
definition. Thus, the fundamental equation (13) is approved: 

{N)p = 2co^ ■ ^, (19) 

where Cq is just a number independent of space dimension, fit giving 1.54 ± 
0.02 in the simulations (Fig.6a). The values of the ratio of c to X/(f) for one-, 
two-dimensional, and mean-field cases are estimated at 3.09, 3.11, and 3.07, 
respectively, which all exhibit the similar value independent of the case of the 
model. The space dimension is irrelevant to the (reduced) prefactor cq. (How 
does one a priori obtain the specific value of the dimensionless quantity Cq 
within the model setting?) 

Thus the above two scaling relations can be unified and reduced to 

W{N; p,p)dN = N-'G (N/{2coXr')pP'') d (N/{2coXr')pP'') , (20) 



with normalization I]*=i^W^(^) = (2coA0^^)pp^^, which corresponding to 
Eq.(17). A scaling function G{x) has a strong drop for a; > 1. To verify this 
law, the simulated data in Figs. 4 and 5 are re-plotted in terms of re-scaled 
coordinates where the x-axis is chosen as x = N/ {2coX(f)~^)pp~^, and the y-axis 
represents y = {2cQX(j)~^) pp~^W (N) . The resulting graph is shown in Fig.6b. 
The data collapse is very good, implying a unified law for school sizes. 

The same numerical results can be obtained with one-dimensional case of the 
mode l . Besi des, the data-collapse in its mean-field case is depicted in Fig. 4 in 
Niwa (|2003i) . Therefore, the unified scaling law, Eq.(20), holds independently 



of the space dimension: space does not influence a linear power-law behavior 
with a crossover to an exponential decay around {N)p. Figure 7 in the next 
section showing the comparison of the empirical school-size distribution with 
simulated distributions in one-, two-dimensional, and mean-field cases directly 
verifies the irrelevance of space in the scaling law for fish school sizes. 



5 The Integrity of Schools 



The cut-off size in the scaling law results from the ability of a school to main- 
tain its integrity over only a certain amount of time: (A^) p varies in proportion 
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Re-scaled population 




0.01 0.1 1 10 

Scaled size 



Fig. 6. Unified scaling law for school sizes, (a) {N) p versus the re-scaled population 
size. The data in Fig.3 are re-plotted in terms of re-scaled coordinate: the abscissa 
represents {X/4>)p/p. Factors X/cp for one- (a), two-dimensional (♦), and mean-field 
(•) cases are estimated at 0.163, 0.469, and 0.985, respectively. The solid line is fit 
with slope 3.08. (b) Simulated school-size distribution plotted on double-logarithmic 
scale (2D). The data in Figs. 4 and 5 are re-plotted with y = cpp~^W{N) as a func- 
tion of the scaling variable x = N/cpp~^, where c = 1.46. The frequency distribu- 
tions are normalized as NW{N) = cp/p. The re-scaled distributions collapse 
onto one another. The data collapse is equivalent to the FSS plo t with A = B = 1 
in Eq. (4) . The solid line is a prediction of the mean- field theory (jNiwal . l2003h . 
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as the coalescence rate, and inversely as the breakup rate [see Eq.(19)]. The 
proposed numerical approach to fish school-size distribution may answer yet 
another challenging question: how schools split. 

Let us analyze the school sizes of Japanese sardine acoustically surveyed in the 
summer in 1982 (the same data set as Fig.l). Transect sampling was carried 
out for a length of ^ = 1211.208 kilometers in a survey area of = 26904 km , 
and 1522 schools (denoted by n) were observed. The average (horizontal) 
interception-length of schools was J = 20.9 meters. We now apply Buffon's nee- 
dle problem to the esti mation of t he dimensional size and the spatial number- 
density of fish schools ( Doi , 1979[ ) . 



Assume the shape of fish schools be horizontally disk-like with the (average) 
diameter R. Let a be the average of inter-school distance, so that the number 
density of schools, 0, is given by a~^. Imagine a disk (instead of a needle) 
dropping on a lined sheet of paper. The probability of the disk hitting one 
of the lines is Pdisk = R/a- Then the average interception- length of schools 
is calculated as £ = n{R/2Y /a. Since we expect to observe one school for a 
transect sample interval of length 6 = ^/fj^, the average length of successive 
detections of schools, the strip of width a and length S contains 1/Pdisk schools, 
i.e. (paS = 1/Pdisk- Accordingly, 6/a = a/R is obtained. Then the average 
diameter of schools and the average inter-school distance are given by 



i?= (/ 1^/5 = 82.6 m, (21) 



and 



3 



7r2 



-M2 = 256.4 m, (22) 

TC 



respectively. Thus the fine-grained spatial number-density of schools, 0, is 



estimated at 15.2 km ^. 



We now return to the coarse-grained model of school aggregation. Let the av- 
erage diameter of schools, R, be the lattice constant. The data are normalized 
to the coarse-grained density of schools, = (pR = 0.104: 

s-'J2m = cj), (23) 



where s = 2^^ sites (~ ). Then the coarse-grained population density is 
estimated at 

p = s-'Y.N^^, = 0.415, (24) 

i 

where unit is half the class width of size histogram (per site). The empirical 
data yield (iV)p = 11.7. 
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Table 4 

Summary of simulations for Japanese sardine schools. Monte-Carlo simula- 
tions have been performed with coarse-grained zones of s = 2^2 sites, $ = 1.74 x 10^ 
individuals, and simulation run = 2^'' time steps. 





P 


<j) ^ 


{N)p - 


half-hfe T ^ 


A 


PT {=P- T) 


ID (A) 


0.0179 


0.105 


11.7 


23.5 


1.15 X 10-2 


0.421 


2D (0) 


0.0516 


0.092 


13.6 


7.55 


4.02 X 10-2 


0.390 


mf (o) 


0.107 


0.109 


10.9 


3.45 (3.29)'= 


9.41 X 10-2 


0.369 



^computed from simulation results at the last of run 
'^Unit is simulation time step. 

'^predicted by (ln2)/(p + 0) with (j) = 0.104 (evaluated value from empirical data) 




1 2 5 10 20 50 100 200 

School size 



Fig. 7. Empirical school-size distribution of sardine observed in the summer in 
1982 (•), same as Fig.l but with normalization Eq.(23), compared with simu- 
lated distributions after 2^'' time steps in one- (a), two-dimensional (0), and 
mean-field (o) cases of the model. The simulated distributions are in good 
agreement with empirical data. Different plots for the empirical observation 
and the numerical simulations fall onto one anoth e r. T he solid line is a pre- 
diction of the mean-field theory [Eq.(ll) in Niwal ()200 3)]. with normalization 



^■[{N)p /Q°°exp[-j;(l -e-^/2)]dx] ^ = 1.31 x 10^ consistent with Eq.(24), where 
{N)p = 11.7. The simulation parameters are summarized in Table 4. Plotted on 
double-logarithmic scale. 

Now we can numerically experiment on the school size distribution with the 
population $ = sp and the breakup rate p listed in Table 4. The breakup 
probability p at each simulation time step is estimated at cp/{N)p, by using 
Eq.(13) with prefactors resulting from simulations (Fig.3). A one-sized school 
of simulated fish {N = 1) is considered as an atomic object, which contain a 
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Fig. 8. Persistency of simulated sardine schools in one- (A), two-dimensional (0), 
and mean-field (o) cases of the model. Shown is a ratio of the number of schools 
maintaining their integrity over a certain amount of time steps of simulation run 
without coalescing or splitting. Broken lines depict the exponential decays with 
rates p+X = 0.0294 (ID), 0.0918 (2D), and 0.201 (mf) fitting the numerical results, 
respectively (unit: per time step). The simulation parameters are summarized in 
Table 4. Plotted on semilogarithmic scale. 

certain number of fish correspond to half the class width of frequency data of 
sardine school sizes observed in the wild. Starting from random initial config- 
urations of eight-sized schools on the lattice, the coarse-grained simulations 
imitate the fine-grained situation of the real- world, as shown in Fig. 7. The re- 
markable consistency between the empirical data and the model's prediction 
unambiguously describes that the space dimension cannot be relevant to the 
scaling exponent of power-law school-size distribution. 

Let us investigate the dynamic properties from the numerical results. There 
are two competing processes in the interacting school system: breakup and 
coalescence. The disintegration behavior of fish schools consists of these two 
processes. The breakup and coalescence rates, therefore, combine into the 
probability of disintegration per unit time per school, which gives the half-life 
T as its reciprocal: 

In 2 

T = -. 25 

p + X ^ ^ 

After time T, the number of schools maintaining their integrity is half of the 
original number (Fig.8). In the mean-field case of the model, the coalescence 
probability per simulation time step is given by the density of schools, A = 0, 
which is tested numerically. The half-life of simulated fish schools is calculated 
in Table 4. The breakup probabihty per half-life, pr {— p • T), obtained for 
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each case (lD/2D/mf) exhibits the similar value independent of the case of the 
model. Therefore , if the half-li f e of sa rdi ne schoo l is experimentally determined 
in the wild [e.g. Lester et al. (198^ or Bavlifi ([l988l) for skipjack tuna], we 
can estimate the breakup rate and the coalescence rate in the wild at pr/T 
and (In 2 — pt)/T, respectively, with numerically predicted value pr- 



6 Discussions 



Data collapse is a way of establishing the scaling in fish school-size distri- 
butions. I have extracted the power-law exponent of size distributions via a 
minimization of a measure to quantify the nature of finite-size scaling collapse, 
in contrast to the 'best-by-eye' data-collapse method. The number W{N) of 
iV-sized schools decays with the truncated power-law form with exponent — 1, 
and the power-law cut-off scales with {N)p, which is a well-defined mean of 
the population distribution among school sizes, NW{N). In order to explain 
the observed scaling property, I have chosen a stationary equilibrium model, 
contrary to a stationary non-equilibrium system of Bonabeau et al.'s (1995, 
1998, 1999) model. It has been found that the scaling law for school sizes of 
pelagic fish is completely independent of the dimension of the space in which 
the fish move: the space dimension is irrelevant to the power-law exponent of 
size distributions. This result is contrary to common knowledge in physics that 
the critical exponents characterizing the scaling behavior observed in phase 
transitions depend on the dimensionality. The model of school aggregation 
does not perfectly mirror the interacting school system in the wild (e.g. the 
size of a school may be relevant to its breakup probability), yet the numeri- 
cally simulated result conform s almost perfec tly to the empirical data. This is 
a consequence of universality ( Stanlev . Il995| ). 



The scaling law supports the view that the power-law distribution of fish 
schools is a self-organized critical phenomenon (Bak, 1996; Jensen, 1998), not 
merely a reflection of an exponential distribution of population among school 
sizes, because only critical processes exhibit data-collapse (Yeomans, 1992; 
Bak et al., 2002, Christensen et al., 2002), known as scaling in critical phe- 
nomena. The interacting school system is naturally attracted to the critical 
value of the spatial number-density of fish schools, without any external ad- 
justment being necessary. In the system there are two competing processes, 
coalescence and breakup, and the critical density depends crucially on the in- 
terplay between the coalescence and breakup time scales. If the school density 
becomes greater (or less) than the critical value, this increased (or decreased) 
density in turn increases (or decreases) the probability of coalescence, leading 
to a shift toward less (or greater) density until it reverts to the critical value. 

The relation between dimensional and social sizes of pelagic fish schools. 
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Eq.(l), gives another scaling law for the scho ol size, which is analogous to 
that used in polymer physics ( de Gennei . 1979I ). What is universal in this law 



is the exponent u: it is the same for all schools, supposed to be independent of 
not only environmental conditions but also species. What is non-universal here 
is the prefactor. It depends on the details of interactions between individuals. 
If we want to understand the properties of schooling configurations, the first 
step is to explain the existence and the value of the exponent v. The second 
step is to account for the constant that multiplies , and this involves deli- 
cate studies on local properties within a school. In the first stage, the school 
size N (in number or biomass) must be measured for different values of di- 
mensional size R and compare them. Finding the exact value of for a given 
school is not an easy task, whereas acoustic surveys for pelagic species are 
extensively performed. 

By contraries, the exponent v may be estimated by hypothesizing FSS in 
the school-size distribution of pelagic species. Assume the scaling relation, 
Eq.(l), and choose the suitable value of v to achieve the best data-collapse 
in the FSS analysis, while a dimensional size histogram for acoustic-survey 
data is transformed into a social size histogram by using Eq.(2). Once the 
scaling relation between dimensional and social sizes of pelagic fish schools is 
established, the precision in the stock assessment will be largely improved. 
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